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چکیده سنگافت یکی از انواع ناپایداری دامنه‌ای است که در اکثر مناطق کوهستانی رخ می‌دهد و موجب صدمات جانی و خسارات سنگین 
به ت سیسات, مناطق مسکونی, جاده‌هاء مزارع و وسایل نقلیه می‌شود. بتابراین برای جلوگیری از صدمات مالی- جانی؛ پیش‌بینی و مدل کردن 
پدیده‌ی سنگ‌افت در دیواره‌های پر شیب برای انجام اقدامات پی شگیرانه ضروری می‌باشد. در اين راستا در مقاله‌ی حاضر به‌منظور ارزیابی 
حطر سنگ‌افت ابتدا بااستفاده از نرم‌افزار [06[01مسیر حرکت سقوط سنگی بلوک‌های جداشده به‌صورت منفرد و جدا از هم در بحرانی‌ترین 
مقطع از ارتفاعات شمال باعتری شهر شاهرود مورد تحلیل قرا رگرفت. سپس به‌منظور درنظ رگرفتن فعل و نفعالات بین ذرات در مسیر 
حرکت, شکل, اندازه و وزن قطعات از روش المان مجزای 3۴۳ استفاده شد. با تحلیل نتایج حاصل از شبیه‌سازی به‌کمک دو نرم‌افزار 
م ذکور, این نتیجه حا صل شد که اگر سنگ‌افت اتفاق بیافتد, تعداد زیادی از قطعات در محدوده‌ی جاده‌ها و مناطق مسکونی توقف دارند که 
این توقف‌ها راه‌بندان در جاده و نعطرات مالی- جانی را در پی حواهد داشت. به‌عبارتی دیگر بانوجه به این‌که بیشترین فاصله‌ی افقی طی‌شده 
ت و سط قطعات ۳٩۳‏ متر و فاصله‌ی اققی مناطق مسکونی ۳۶۰ متری می‌با شد, لذا احتمال حادثه در اثر سنگ‌افت بسیار بالا است. در ادامه‌ی 
پیشنهاداتی از جمله نصب توری‌های فلزی» ضربهگیرها, کاهش انرژی جنبشی در قطعات, تغییر در ابعاد قطعات, ایجاد تغییر در هندسه‌ی 


شیب جهت کاهش خعطر و حسارات ناشی از سنگ‌افت در مقطع موردنظر ارائه شده‌است. 
واژه‌های کلیدی سنگافت. روش توده‌ای- کلوخه‌ای» نرم‌افزار 1800111 روش المان مجزاء نرم‌افزار ۳۳0۴ 
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مقد مه 

تحلیل پایداری بلوک‌های سنگی منفصل واقع بر شیب‌ها 
و دامنه‌ها که تحت عنوان سنگ‌افت نامیده می‌شود امروزه 
سبب شده‌است تا بسیاری از مهندسان و طراحان به 
پیش‌بینی رفتار این پدیده در شیب‌های سنگی روی 
آورند. عوامل متعددی در رخداد سنگ‌افت نقش دارند 
که از جمله‌ی آنها می‌توان به افزایش آب حفره‌ای در 
وهی نگ فرهایشن ممانم خر که هر زو کی در 
سنگ. زمین‌لرزه» اثر گوه‌ای ریشه‌ی درختان» تجزیه‌ی 
شیمیایی سنگ. فرآیند سریع یخ‌زدن و آب‌شدن اشاره 
کرد. اثرات ناپایداری دامنه‌ها به‌صورت حرکات توده‌ای 
در مقیاس کوچک تا وسیعی می‌باشد. حرکات دامنه‌ای 
ممکن است جزئی و منحصر به ریزش یک قطعه سنگ 
منفرد بوده یا اين‌که بسیار بزرگ و فاجعه‌آفرین باشد. 
مکانیسم حرکتی رفتار سنگ بعد از جدا شدن از سطح 
شیب‌دان در حقیقت ترکیبی از چهار فاز مختلف و اصلی 
لغزش. سقوط آزاد. چرخش و جهش می‌باشد. بنابراین 
می‌توان با تعیین درست مسیر ریزش قطعات که به‌صورت 
ارتفاع و طول مسیر. سرعت و انرژی قطعات طی ریزش 
بیان می‌شود. برای تعبین محل نصب سازه‌های حفاظتی 
مناسب اقدام کرد. انتهای مسیر حرکت قطعات به عواملی 
نظیر زبری» شیب. ارتفاع ریزش مربوط به دامنه و سه 
فاکتور استحکام شکل و حجم قطعات بستگی دارد. در 
طول سالیان اخیر در ارتباط با پدیده‌ی سنگ‌افت 
پژوهش‌های متعددی انجام شده که در ادامه به مهم‌ترین 
آنها اشاوع تقو اهد: سل 

اززونی ([۸2200) و همکاران بااستفاده از مدل 
ریاضی و براساس مکانیزم جسم صلب و آنالیز آماری در 
یک فضای دوبعدی به بررسی بعضی از پارامترهای 
کاربردی سنگ‌افت پرداختند [1]. آززونی و فریتاس 


( و1)2ز۳ ۸220) بااستفاده از ازمایش‌های برجا با 
پرتاب قطعات سنگی و فیلم‌برداری از تمام مراحل 
حرکت آنها؛ پارامترهای مناسبی را برای تحلیل‌های 
ریزش سنگی ارزیابی کردند. اين پارامترها شامل ضریب 
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و غلتشی و ارتفاع جهش می‌باشد. آنها هم‌چنین تأثیر 
شکل و ابعاد بلوک‌ها را روی سرعت غلتشی و اثر وجود 
گودال را برروی مسیرهای ریزش بررسی کردند [2]. 
نیکت (600:) و همکاران بااستفاده از روش‌های تحلیل 
و آزمایشگاهی» تأثیر شبکه‌های فلزی (-9اک-ناه۸ 
۷1706 در توقف سنگ‌افت را بررسی کردند. آنها نشان 
دادند که این شبکه‌هاء دربرابر نیروی دینامیکی افت سنگ 
انعطاف‌پذیری بالایی دارند [3]. کوزتی (8ع0072) و 
همکاران به‌منظور بررسی سقوط سنگ در مقیاس 
ناحیه‌ای و منطقه‌ای و به‌صورت سه‌بعدی برنامه‌ی 
0۲ را طراحی کردند که قابلیت مدل کردن 
سنگ‌افت را به‌صورت سه فاز سقوط آزاده غلتش و 
جهش دارد. داده‌های ورودی این برنامه در محبط 
سیستم اطلاعات جغرافیایی (018) تهیه شد که شامل 
ژئومورفولوژی. سرعت ابتدایی. ضرایب دینامیکی 
ارتجاعی و اصطکاک می‌باشند. برنامه‌ی مذکور دارای 
قابلیت مدل کردن منبع سقوط سنگ به سه حالت خطی» 
نقطه‌ای و ناحیه‌ای می‌باشد که نتایج حاصل به دو شکل 
دوبعدی و سهبعدی در محبط 015 ارائه می‌شوند [14. 
دوارت و مارکوئیز (عتعدظ مها تانشز 
فاکتورهای محیطی و لیتولوژیکی را بر ایجاد ریزش‌های 
مننگی در مقیاس ناحیه‌ای و بااستفاده از نهشته‌های 
واریزه‌ای در محبط 015 بررسی نمودند. در این مطالعه 
از واریزه‌های جدید که توسط پوشش گیاهی دربرگرفته 
نشده‌اند. به‌عنوان نشانی برای فعالیت سقوط سنگ 
استفاده شده‌است. واریزه‌هایی که توسط گیاهان پوشیده 
شده‌اند. دال بر پایداری شیب‌ها و عدم فعالیت آنها 
می‌باشد. در تحقیق مذکور نقشه‌ی حوضه‌های سقوط 
سنگ با نقشه‌های متغیرهای محیطی (دما و باران), 
ساختمانی و ژئومورفولوژی هم‌پوشانی شده‌است. نتایج 
حاصل حاکی از آن است که متغیرهای محیطی تأثیرات 


عمده‌ای بر فعالیت دامنه‌ها دارند. به‌حصوص که این 
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مرها در یرلوفی‌های: متتاف: تایرات تفای ار 
ایجاد می‌کنند. نتایج به‌دست‌آمده فعالیت شدید سقوط 
سنگ در نقاط مرتفع‌تر را تأیید می‌کند. در مناطقی که 
لیتولوژی‌های مقاوم و نامقاوم به‌طور متناوب قرار دارند. 
سقوط سنگ‌ها بیشتر مشاهده می‌شود [5]. چائو (دهه) 
و همکاران بااستفاده از یک سری آزمایش‌های 
آزمایشگاهی ضرایب واجهشی برای بلوک‌های 
کروی‌شکل در شیب‌های سنگی را به‌دست آوردند. آنها 
هم‌چنین دریافتند که ضرایب واجهشی با افزایش زاویه‌ی 
شیب افزايش می‌یابد. هم‌چنین به‌طور کلی نسبت انرژی 
سیستماتیک چرخشی به انتقالی با افزایش زاویه‌ی شیب 
تا ۰ درجه افزايش و سپس تا زاویه‌ی شیب ۷۰ درجه 
کاهش می‌یابد [6]. بیلی‌فرد (۵ننللنه8) و همکاران 
نقشه‌ی خطر سقوط سنگ را بااستفاده از عوامل 
توپوگرافی. ژئومورفولوژی و زمین‌شناسی و بااستفاده از 
5 تهیه نمودند [7]. کوزتی و همکاران یک روش 
براساس افزایش حوادث ناشی از سنگ‌افت بااستفاده از 
تاریخچه و حجم سنگافت و برپایه‌ی نتایج فیزیکی و 
مدل فضایی شبیه‌سازی سنگافت برای تعیین خطر 
سنگ‌افت ارائه کردند. آنها اطلاعات مربوط به محل و 
نوع اقدامات مقابله‌ای دربرابر سنگ‌افت را بیان و این 
روش را در دره‌ی رودخانه نرا در منطقه‌ی 0712 
ایتالیای مرکزی مورد استفاده قرار دادند [8]. چائو و 
همکاران بااستفاده از نرم‌افزار 015 به پهنه‌بندی خطر 
سقوط سنگ در هنگ‌کنگ پرداختند. آنها در تحقیق خود 
از چهار لایه‌ی اطلاعاتی زمین‌شناسی. مقدار و جهت 
شیب و طبقات ارتفاعی به‌عنوان عوامل اصلی در وقوع 
پدیده‌ی سقوط سنگ در این منطقه استفاده کردند [9]. 
پرت (۳۵۳۵) و همکاران به مطالعه‌ی فعالیت مکانی و 
زمانی ریزش در مناطق جنگلی کوه‌های آلپ سوئیس 
پرداختند و نتیجه گرفتند وقوع ریزش فصلی در طول 
خواب درختان. در اوایل بهار به اوج می‌رسد. علاوه‌بر 
این» وقوع ریزش با درجه‌ی حرارت تابستان و زمستان 
در ارتباط است. ولی هیچ هم‌بستگی با بارش سالانه یا 
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فصلی ندارد [10]. سابویا (و5200) و همکاران برای 
ارزیابی ناپایداری‌های دامنه در ریودوژانیروی برزیل. از 
منطق فازی استفاده کردند. آنها به اين نتیجه رسیدند که 
در این مدل چون به‌صورت صفر و یکی نیست. نتایج 
بهتری نسبت به سایر مدل‌ها به‌دست می‌آید [11]. 
باترسون (826۲500) و همکاران در مطالعه‌ای ارزیابی 
خطرات زمین‌شناسی از جمله سنگ ریزش را در دره‌ی 
هامبر در شمال شرق آوالون در کانادا انجام دادند. و به 
اين نتیجه دست یافتند که در شیب‌های ۲۲ تا ۳۰ درجه 
بیشترین ریزش‌های منطقه اتفاق افتاده‌است [12]. توپال 
(1001) و همکاران به بررسی سنگ‌افت در شهر افیون 
ترکیه پرداختند. قلعه‌ی افیون یک مکان توریستی در شهر 
افیون ترکیه است که در یک تبه‌ی شیب‌دار قرار 
گرفته‌است و در مجاورت قلعه سکونتگاه‌هایی وجود دارد 
که در معرض خطر سنگ‌افت قرار دارند. آنها بااستفاده از 
نرم‌افزار 1800111 به بررسی سنگ‌افت و محاسبه‌ی ارتفاع 
جهش. انرژی جنبشی. سرعت سنگ‌ها و نقطه‌ی توقف 
بلوک‌ها در طول مسیر پرداختند. نتایج تجزیه و تحلیل‌ها 
نشان داد که برای مناطقی که مستعد خطر ریزش سنگین 
هستند. اقداماتی نظیر پیچ‌سنگ و نرده‌های محافظ قابل 
پيشنهاد است [13]. لان (عمآ) و همکاران برای تسهیل 
ارزیابی خطرات ریزش سنگی در محیط 15 
زیربرنامه‌ای را تحت عنوان ۸7٩۲‏ 00191 ارائه 
نمودند که ترکیبی از یک مدل‌سازی سه‌بعدی فرآیند 
وی مکی با مد اقا 6019 اس 12 کوش 
(نمنصمعدن) و همکاران با انجام ۲۰ آزمایش سقوط 
سنگ, دید تازه‌ای در مورد قطعه‌قطعه شدن بلوک‌های 
سقوط کرده به‌خصوص بلوک‌های دارای فولیاسیون 
هنگام برخورد با سطح زمین ایجاد کردند. اين آزمایش‌ها 
در معدنی در شمال‌غرب ایتالیا برروی دو نوع سنگ با 
شیستوزیته‌ی متفاوت صورت گرفت [15]. چوی (ن0)) 
و همکاران بااستفاده از مدل بارتون و نرم‌افزار 100011 
به حل مشکل خطر ریزش ایالت سوک گای‌اوراک در 
کره‌ی جنوبی پرداختند و بیان کردند ۸۰ درصد 
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تخته‌سنگ‌های در حال سقوط دارای وزن حدود ۱۰۰ تا 
۰ کیلوگرم هستند. میانگین سرعت آنها ۸/۶ متر در 
ثانیه است و حداکثر انرژی جنبشی این سنگ‌ها ۱۲/۵ 
کیلوژول است [116. تونسول اوغلو و زورلو 
( بلعمتاونمت؟ بلبم2 ) به‌ کمک مشاهدات میدانی و 
بااستفاده از نرم‌افزار 111 نقشه‌ی پهنه‌بندی وقوع 
ریزش قلعه‌ی اورتا حصار در منطقه‌ی کاپادوکیه در کشور 
ترکیه را ترسیم کردند [17]. آندریس و لکاس ( «0نهنهه 
14) در تحقیق خود به بررسی و آنالیز پدیده‌ی 
سقوط سنگ در قسمتی از جزیره‌ی ۹۵۵10201 در کشور 
یونان پرداختند. در این بخش از جزیره» سقوط سنگ در 
اطراف جاده خطراتی را برای وسایل نقلیه و ساکنان 
محلی ایجاد کرده بود. آنها بااستفاده از نرم‌افزارهای 
۷ و 50011 منطقه‌ی مورد مطالعه را پهنه‌بندی 
و آن را به سه زون با خطر کم متوسط و زیاد تقسیم 
کردند [18]. رممر (عصصمع) و همکاران در رشته 
کوه‌های آلپ از مدار اتریش تا سوئیس بااستفاده از مدل 
11 به بررسی خطر ریزش پرداختند [19]. هوانگ 
(عمهت۳۳) و همکاران با مطالعه درمورد اثرات بارندگی‌ها 
بر گسیختگی‌های کم‌عمق دامنه‌ای به‌روش تجربی بدین 
نتیجه رسیدند که با پیدايش هرگونه گسیختگی در پای 
دامنه تحت هر عاملی. گسیختگی‌های متوالی به بخش 
بالادست دامنه. مترایتت هی کند .ی ارقاط. تدحانتکی آبین 
عوامل درونی و بیرونی در شکل‌گیری گسیختگی و 
ریزش سنگ‌ها در جاده‌ها وجود دارد [20]. ولکوین 
(۷۵۳۵8) و همکاران در مطالعه‌ای به بررسی کلی. 
انواع مدل‌های مطرح در پهنه‌بندی ریزش و نیز اقدامات 
حفاظتی درراستای مقابله با این خطر پرداختند [21] 
اسپاداری (ته0ه50) و همکاران؛ ارزیابی آماری از 
محدوده‌ی انرژی سنگافت برای پنج شرط زمین‌شناسی 
ملک قر استولا زا بروبین کردنکه ناما لین مار 
تصادفی دریافتند که برای دو مطالعه‌ی موردی سنگ‌افت. 
میانگین انرژی ضربه ۲۰۰ تا ۳۰۰ کیلوژول و برای سه 
مورد دیگر کمتر از ۱۰۰ کیلوژول می‌باشد که ۹۵ درصد 
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تحلیل بدیده‌ی سنگ‌افت در دیواره‌ی پرشیب بااستفاده از مدل‌سازی... 


این مقدار انرژی باید برای طراحی سازه‌های حفاظتی 
دربرابر سنگ‌افت درنظر گرفته شود [22]. یوسف 
(/۷0۵:0) و همکاران در تحقیق خود به بررسی و آنالیز 
سنگ‌افت در کوه النور در عربستان پرداختند. کوه النور 
یکی از مناطقی است که به‌دلیل وجود غار هیراء پیست 
اسکارپمنت و مناطق شهری اطراف آن به‌علت ریزش 
شیک سیان ترشطر استت, میفان ادست. تافتنل که غیت 
تند. هوازدگی. فرسایش, اثرات انسانی و بارندگی از 
دلایل سنگ‌افت در اين منطقه است. بنابراین آنها سعی 
کردند که بااستفاده از سنجش از راه دور نرم‌افزارهای 
عددی و مطالعات میدانی» مناطق ناپایدار را شناسایی 
کنند. تجزیه و تحلیل نشان داد که ریزش سنگ در 
منطقه‌ی ۱ و ۲ کوه النور ممکن است به مناطق شهری 
برسد» درحالی که ریزش سنگ در منطقه‌ی ۳ به منطقه‌ی 
شهری نخواهد رسید و ریزش سنگ در امتداد جاده‌ی 
پیست آلنور تأثیر زیادی برروی گردشگران دارند و 
بایستی در هر ۳ منطقه اقدامات پیشگیرانه‌ای انجام شود 
[23]. کسکین (90ع) و همکاران به بررسی خطر 
سنگ‌افت در جاده‌ای می‌پردازند که در پایین‌دست آن 
مناطق مسکونی عبور می‌کند. تجزیه و تحلیل سنگافت 
در جاده‌ی زونگولداک- کلیملی در ترکیه براساس 
تحقیقات میدانی» مطالعات آزمایشگاهی و بااستفاده از 
نرم‌افزار 180681 انجام شده‌است. آنها به این نتیجه 
رسیدند برای جلوگیری از خطرات ساختمان‌های 
مسکونی در پایین‌دست. از آنجا که تغییر هندسه‌ی شیب 
امکان‌پذیر نیست. سد بتنی در کنار جاده را به‌عنوان یک 
اقدام پیشگیرانه توصیه کردند [24]. هی (۲16) و همکاران 
با بررسی‌های دقیق میدانی. نقشه‌برداری زمین‌شناسی» 
تجزیه و تحلیل حرکتی و مدل‌سازی عددی خطر ریزش 
سنگ‌ها را در شهر میانسی از شهرستان ونچوان, استان 
سیچوان چین تحت تأثیر زلزله بررسی کردند. آنها به این 
نتیجه رسیدند که سنگ‌های درزه‌دار پس از زلزله سبب 
جدایش شده و علاوه‌بر ریزش در جاده سبب تحریک 


سنگ‌های بدون درزه در پایین‌دست می‌شوند. آنها 
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هم‌چنین اثر بارندگی برروی سنگ‌افت را مورد بررسی 
قرار دادند [125. 

مطابق با تحقیقات انجام‌شده تاکنون مطالعات 
زیادی در این خصوص انجام شده‌است. اما از آنجا که 
نرم‌افزار 13008011 تنها مسیر حرکت را شبیه‌سازی می‌کند 
و عواملی نظیر شکل, اندازه. وزن قطعات و هم‌چنین فعل 
و انفعالات بین ذرات در مسیر حرکت (ذرات به‌صورت 
منفرد و جدا از هم) را درنظر نمی‌گیرد و شبیه‌سازی 
نمی‌کند از واقعیت سنگ‌افت فاصله دارد. درصورتی که 
در واقعیت. قطعات در حال ریزش و سقوط برروی 
یکدیگر تأثیرگذار هستند. بنابراین در اين تحقیق به‌منظور 
تطابق بیشتر با محیط واقعی منطقه‌ی شاهرود. و برای 
پوشش و رفع نواقص نرم‌افزار الل0 از نرم‌افزار 
پرقذوت: ۳۳0 انشفاده شدمانست تا علاومر شناسایی 
سرعت و مسیر سنگ‌هاء. فعل و انفعالات بین ذرات و 
دیواره‌ها و بین ذرات با یکدیگر در حین حرکت در 
پایین‌دست درنظر گرفته شود تا به شرایط موجود در 
طبیعت نزدیک‌تر باشند و بتوان نتایج قابل قبول و 
مطمئن تری را از شبیه‌سازی به‌دست آورد. بنابراین در اين 
مقاله بااستفاده از دو نرم‌افزار مذ کور به‌منظور کامل تر شدن 
اطلاعات در جهت بررسی خطر سنگ‌افت. به شبیه‌سازی 
مکانیسم ریزش سنگی. تعیین میزان جابه‌جایی توده‌های 
ریزشی و هم‌چنین بررسی راهکارهای مقابله و جلوگیری 
از این ریزش‌ها و خطرات ناشی از آن‌ها در منطقه‌ی 


شاهرود پرداخته شده‌است. 


معرفی نرم‌افزار [1800121 
نرم‌افزار 100/011 یک برنامه‌ی تحلیلی- آماری دوبعدی 
ست که توسط استیونز (516۷605) برای شبیه‌سازی سقوط 
تیگ ایجاد شده‌است. استفاده از این روش تحلیل 
حتمالی و متغیرهای تصادفی در متغیر جرم شرایط اولیه 


و سرعت سنگ و هم‌چنین در متغیرهای ضرایب 


رتجاعی. زاویه‌ی اصطکاک و موقعیت رئوس شیب‌ها 
در شبیه‌سازی این برنامه نتایج آماری مفیدی در اختیار 
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قرار می‌دهد و مشکلات مرتبط با عدم قطعیت داده‌های 
ورودی را کاهش می‌دهد. خروجی‌های این برنامه 
به‌صورت منحنی‌ها و نمودارهای توزیع ماکزیمم سرعت.؛ 
انرژی جنبشی, ارتفاع پرش و مکان توقف سنگ در 
سرتاسر طول شیب می‌باشد [26]. از نتایج این نرم‌افزار 
می‌توان برای طراحی دیوارهای حائل استفاده کرد. 
تغییرات جزئی در مقدار ضریب ارتجاعی تأثیرات کاملا 
آشکاری در مقدار متغیرهای خروجی این برنامه از جمله 
ارتفاع پرش سنگ‌ها و مسافت طی‌شده دارد و محاسبه‌ی 
مقدار دقیق آن نتایج شبیه‌سازی را صحیح‌تر می‌کند که 
براساس آن می‌توان محل, اندازه و استحکام مورد نیاز 
دیوار حائل را تخمین زد [26] 

این برنامه از سه الگوریتم برای آنالیز و محاسبه‌ی 
حرکت سنگ استفاده می‌کند که عبارتند از الگوریتم ذره 
الگوریتم جهش و الگوریتم لغزش. الگوریتم ذره همه‌ی 
پارامترهای شبیه‌سازی را دریافت می‌کند و شرایط اولیه 
را برای الگوریتم سقوط و الگوریتم پرتاب فراهم 
می‌نماید. الگوریتم پرتاب برای محاسبه‌ی حرکت سنگ 
در حالت جهش و الگوریتم سقوط برای محاسبه‌ی 
حرکت در زمانی که سنگ با سطح شیب‌دار برخورد و 
به‌سمت پایین حرکت می‌کند. استفاده می‌شود. از تلاقی 
یی با کل ی یرومیت که ور 
لحظه‌ی برخورد و حل معادله‌ی درجه‌ی دوم می‌توان 
سرعت. انرژی و ارتفاع پرش را در هر نقطه تعیین نمود 
[26]. 


معرفی نرم‌افزار المان مجزای ۳1 
نرافزار ۳6 از مجموعه نرم‌افزارهای گروه 17۸50۸ 
که اشنم سای نی تقد فا نی زب 
نرم‌افزار از روش المان مجزا برای شبیه‌سازی مدل‌های 
خود استفاده می‌کند. در این روش توده‌سنگ به‌صورت 
یک محیط ناپیوسته و برمبنای روش اجزای گسسته و 
ذره‌ای درنظر گرفته می‌شود که رفتار مواد را بااستفاده از 


محبط ذره‌ای و دایره‌ای در حالت دوبعدی و سه‌بعدی 
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تحلیل می‌کند. بنابراین تأثیر عوارض ساختاری نظیر 
درزه‌های سنگی, لایه‌بندی. گسل‌ها و انواع ناپیوستگی‌ها 
قابل بررسی است. این روش باتوجه به ماهیت ناپیوسته 
بودن محبط‌های سنگی روش کاملا مناسبی برای بررسی 
تغییرشکل و شکست توده‌های سننگی می‌باشد. در این 
روش برای شبیه‌سازی مسیرهای توده سنگ و محاسبه‌ی 
انرژی سنگ در حال سقوط و روابط برخورد بین عناصر 
استفاده می‌شود. نرم‌افزار ۳۳6 حرکت و تعامل با سطح 
مقطع دایره‌ای (219) و يا کروی (31) ذرات و عناصر 
دیوار را بااستفاده از قوانین حرکت و نیروی جابه‌جایی 
مدل‌سازی می‌کند. در طول محاسبه. تماس بین ذرات با 
ذرات يا ذرات با دیوارها به‌طور خودکار شناسایی 
می‌شود. ذرات ممکن است در نقاط تماس باهم پیوند 
داشته باشند و این پیوند به‌سبب برخورد می‌تواند شکسته 
شود. بنابراین می‌توان بااستفاده از دستورات نرم‌افزار 
انواع حرکت (جهشی: لغزشی» غلتشی و سقوط 
آزاد یک توده سنگ منفرد) و تعامل بین بلوک‌ها را در 
1 


موقعیت جغرافیایی و زمین‌شناسی محدوده‌ی 
مطالعاتی 
شهر شاهرود بزرگ‌ترین شهرستان استان سمنان می‌باشد 
که از شمال به ارتفاعات سلسله جبال البرز و شهرستان 
گنبد کاووس و گرگان» از شرق به شهرستان سبزوان از 
جنوب به دشت کویر و استان اصفهان و از غرب به 
شهرستان دامغان مرتبط می‌شود. شاهرود به‌لحاظ موقعیت 
ممتاز و ویژه‌ی جغرافیایی و قرار گرفتن بین رشته‌کوه‌های 
البرز و سرزمین‌های کویری همواره از جمله جاذبه‌های 
گردشگری ایران بوده‌است. از آنجایی که اين منطقه یک 
مکان تفریحی در شمال شهر شاهرود می‌باشد هر ساله 
گردشگرانی برای بازدید از اين منطقه به این شهر سفر 
می‌کنند اما از آنجا که این محدوده در معرض برخورد 
قطعات سنگی است در مواردی ریزش توده‌سنگ در 
جاده‌ها سلامتی گردشگران را به خطر می‌اندازد. هم‌چنین 
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مزار شهدای گمنام در پایین‌دست منطقه. فضای سبز واقع 
در شمال هتل جهانگردی. شماری از منازل مسکونی 
پایین‌دست. بخشی از فضای سبز حاشیه‌ی جنوبی محل 
تفریحی آبشار و حاشیه‌ی باختری جاده‌ی کاج نیز در این 
محدوده قرار دارد؛ بنابراین بررسی پدیده‌ی سنگ‌افت و 
راهکارهای مقابله با آن برای این منطقه بسیار حائز 
اهمیت می‌باشد. 

منطقه‌ی مورد مطالعه ارتفاعات شمال باختری شهر 
شاهرود است. برای مطالعه‌ی دقیق‌تر نقشه‌ی زمین‌شناسی 
بزرگ‌مقیاس منطقه (۱/۱۰۰۰۰) از عکس‌های هوایی 
انکفاده تقیله‌اس:, کل 1(۰ فان وشگتاس 
۰ منطقه‌ی مورد مطالعه را نشان می‌دهد. 

قله‌ها و بخش مرتفع ارتفاعات. متعلق به سازند 
آهکی صخره ساز لار است. در زیر آهک لار. سازند 
یجای: قارب آشکن عارین او مارن): قران داز که 
نسبت‌به فرسایش مقاومت کمتری دارد. رخنمون کوچکی 
از شیل زغال‌دار شمشک نیز در زیر سازند لار به‌صورت 
فرسایش یافت می‌شود. در همبری سازندها اختلاف 
فرسایش سبب شده‌است که آهک‌های لار به‌صورت 
معلق درآیند. علاوه‌بر عامل فرسایش, حضور دو گسل 
موازی با روند ارتفاعات یکی در مرز سازند لار با دلیچای 
و شمشک و یکی درون سازند لار سبب به‌وجود آمدن 
پرتگاه‌هایی به‌ارتفاع ۳ تا ۶ متر شده‌است. تکتونیزه بودن 
سازندهاء هوازدگی ناشی از اختلاف دمای شبانه‌روزی و 
وجود شمار زیاد روزهای یخبندان در طول سال و فرایند 
ذوب و انجماد آب در درزه و شکاف‌ها سبب شده‌است 
تا پرتگاه‌های موجود استعداد بالایی برای ریزش قطعات 
شتگی دافقة ماففد: کته (۲) نان از واریرهه رانر 
محل نشان می‌دهد. در برخی از مناطق نیز قطعاتی بسیار 
بزرگ به ابعاد چند ده مترمکعب دیده می‌شود. 


سال سی و چهارم» شمارة یک» ۱۶۰۰ 


هادی فتاحی- میلاد صادق یکواری- حسین قاثدی ۱۱ 


شکل ۳۲ واریزه‌های موجود در دامنه‌ی مشرف به جاده و مناطق مسکونی 


همان‌طور که ذکر شد منطقه‌ی مورد مطالعه متعلق جمله رخداد کوه‌زایی. چین‌خوردگی و گساش‌های 
به ایالت زمین ساختی البرز شرقی است. البرز شرقی در متعدد بوده‌است. در ورقه‌ی زمین‌شناسی ۱/۱۰۰۰۰۰ 
طی ادوار زمین‌شناسی متحمل رویدادهای مختلفی از شاهرود چین‌خوردگی‌ها شامل ساختمان‌های تاقدیسی و 


سال سی و چهارم» شمارة یک ۱۶۰۰ نشریه مهندسی عمران فردوسی 


ناودیسی با روند شمال شرق- جنوب غربی و يا شرقی - 
غربی و از انواع متقارن یا نامتقارن می‌باشد. سیستم اصلی 
کی ها هضور نگ ها و شیس کسا هاش مداد لش 
است. راندگی‌های موجود دارای امتداد شرقی- غربی 
است و حرکات قابل توجهی را از شمال به جنوب انجام 
داده‌اند. این راندگی‌ها در جاپی که شیب گس به صفر 
نزدیک شده‌است. سفره‌های رورانده را ب‌وجود آورده‌اند. 
در محدوده‌ی مورد مطالعه نیز واحدهای مختلف 
چینه‌شناسی دارای چین خوردگی‌های واضح و مشهودی 
می‌باشد. از آنجا که توده‌سنگ به‌صورت یکپارچه نیست 
و دارای درزه و شکستگی‌های متعدد می‌باشد و از طرفی 
در تنش‌های پایین که در اعماق کم وجود دارد. رفتار 
توده سنگ توسط شکستگی‌ها کنترل می‌شود. لذا شناخت 
خحصوصیات درزه‌ها در ارزیابی پایداری شیب‌های سنگی 
ضروری می‌باشد. به این منظور. ارتفاعات شاهرود که 
دارای درزه‌های مستعد ریزشی می‌باشد با رس و سیلت 


پر شده‌اند. 


نتایج حاصل از مدل‌سازی با روش توده‌ای- 
کلوخه‌ای 100۲211 
به‌منظور بررسی وضعیت پایداری دامنه در منطقه‌ی مورد 


مطالعه باتوجه به نقشه‌ی توپوگرافی و بازدیدهای 


تحلیل بدیده‌ی سنگ‌افت در دیواره‌ی پرشیب بااستفاده از مدل‌سازی... 


ی ایب فان رای ای فصو تفا 
منشاً ریزش سنگی, اصلی‌ترین معیار برای بحرانی بودن 
شرا در که ند ورگ هی ی رسای 
شیب تند و گاهی معلق یا شیب منفی. حضور قطعات 
ریزشی بسیار بزرگ در پای دامنه پا برروی دامنه و وجود 
مناطق مسکونی در فاصله‌ی نزدیک به دامنه می‌باشد. با 
این اوصاف مقطعی مورد ارزیابی قرار گرفت که به‌دلیل 
وجود منطقه‌ی مسکونی در پایین‌دست این مقطع» در 
صورت ریزش سنگ‌افت. شرایط بحرانی را برای ساکنان 
آن ایجاد می‌کرد. خصوصیات درزه‌ها در مقطع مدنظر که 
سبب ریزش سنگافت شده مطابق جدول (۱) و (۲) 
تفه 

همان‌طور که قبلا اشاره شد. از نرم‌افزار 100۳211 
برای شبیه‌سازی پدیده‌ی سنگ‌افت در مقطع مورد نظر 
استفاده شده‌است. این مقطع باتوجه به اندازه‌گیری‌های 


انجام شده به‌صورت شکل ۳ می‌باشد. 


جدول ۱ خحصوصیات درزه‌های اصلی در مقطع مورد نظر 


دستق. | مبانکنم هت فاصله‌داری (60) 2 حداکثر 
درزه | شیب | شیب | ماکزیمم | میانگین | مینیمم ِِ بازشد گی(00) 
۱ ۶ | ۱۳۷ | ۱۸۰ ۹۶ ۳۰ ۳ 1 
۲ 11 ۳0۰ ۹ ۷۸ ۳ رس ۱ 
۳ ۸ ۳۱ ۱۳۷ 5:1 ۱1 سیلت ۱۰ 
3 ۷ ۳:۲ ۱۸۳ ۷ ۷/۳ رس ۲ 
۵ ۸ ۳۲ ۱۸۰ ۹ ۳۰ رس ۳ 


جدول ۲ خصوصیات مقاومت برشی درزه‌های اصلی در مقطع مورد نظر 
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سال سی و چهارم. شمارة یک» ۱۶۰۰ 


هادی فتاحی- میلاد صادق یگواری- حسین قائدی 


ضریب زبری (1156) 
دسته درزه « 
ماکزیمم | میانگین | مینیمم 
۱ ۳۱ ۸ ۳۳ 3 
1 ۳۱ ۹ ۷ ۵ 
۳ ۳۱ ۸ ۷ ۵ 
3 ۳۱ ۸ 23 3 
۵ ۳۱ ۸ ۷/۲ 1 
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شکل ۳ نیم‌رخ مقطع مورد نظر (واحد محورهای مختصات به متر است) 


در شکل (۳) نقاط منشأً ریزش تعیین شده‌اند و برای 
بخش‌های مختلف دامنه مطابق جدول (۳). ضرایب 
ارتجاعی. ضریب اصطکاک و زبری متناسب با 
ویژگی‌های آن تعیین و وزن بزرگ‌ترین قطعات در هر 
مقطع بااستفاده از برداشت‌های صحرایی برپایه‌ی 
فاضهداز ی دروه‌ها مقتصی شاه‌اسه: 

مطابق شکل (۶) در اين نرم‌افزار برای شبیه‌سازی 
۰ ریزش در طی هر بار اجرای برنامه تنظیم شده‌است 
و بزرگ‌ترین قطعه در مقطع دارای وزن معادل ۳۸۰۰ 
کیلوگرم می‌باشد. 

باتوجه به تعداد سنگ‌های در طول مسیر سقوط 
می‌توان مشاهده کرد که پس از ویزش محل #وقف تعداد 
زیادی از سنگ‌ها در فاصله‌ی ۲۸۰ متری تا ۳۶۰ متری که 


محل عبور جاده است می‌باشد. هم‌چنین تعدادی از 


سال سی و چهارم شمارهُ یک, ۱۶۰۰ 


قطعات در مسیر جاده توقف دارند. درصورت سنگ‌افت. 
این توقف‌هاء باعت بسته شدن جاده و ایجاد حوادث و 
قآ ان اند بر ماکان انم مقاظی هی تیه 
همان‌طور که در شکل (۵) مشخص است. تغییرات در 
سرعت اولیه و وزن قطعات تأثیر قابل توجهی در تعداد 
و محل توقف قطعات ندارد. 

در این منطقه و مطابق شکل (1) باتوجه به ارتفاع 
پرش قطعات سنگی در طول مسیر می‌توان محدوده‌هایی 
که قطعات لغزش, غلتش با جهش دارند را مشخص کرد؛ 
از نمودارهای به‌دست‌آمده به‌کمک نرم‌افزار مشخص 
می‌شود که قطعات در حد فاصل ۹۵ تا ۱۱۰ متری با 
جهشی حدود ۰/۷۵ متر و در فاصله‌ی ۲۵۰ متری جهشی 
بالای ۷ متر دارد. ابتدا و انتهای مسیر حرکت: قطعات 
به‌صورت لغزش يا غلتش حرکت می‌کنند. هم‌چنین 
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۱۰۴ تحلیل بدیده‌ی سنگ‌افت در دیواره‌ی پرشیب بااستفاده از مدل‌سازی... 


باتوجه به تغییرات در وزدن و سرعت اولیه و نمودارهای 
حاصل از آنها می‌توان به این نتیجه رسید که اين تغییرات 


جدول ۳ ویژگی‌های مواد سازنده‌ی دامنه‌ی مورد نیاز در ورودی برنامه 


بخش بالایی دامنه بخش پایینی دامنه 
۳ ۷ سنگ آهک مارنی پوشیده با مقدار زیادی 
نوع سنگ‌های دامنه سنگ اهک‌های توده‌ای و مقاوم لار 2 ا 
سنگ‌ریزه و کمی پوشسس گیاهی 
زبری دامنه (انحراف معیار از شیب میانگین) ۲ ۲ 
زاویه‌ی اصطکاک دامنه (درجه) ۳۰ ۳۰ 
ضریب ارتجاعی عادی (بدون واحد) ۰/۵ ۳۳ 
ضریب ارتجاعی مماسی بدون واحد) ۰/۹۵ ۰/۸۷ 
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150017 
14501 
14001 
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شکل ۶ مسیرهای طی‌شده‌ی قطعات در مقطع در طی ۰ بار شبیه‌سازی 
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شکل ۵ مسافت طی‌شده 


نشریه مهندسی عمران فردوسی سال سی و چهارم» شمارة یک, ۱۶۰۰ 


هادی فتاحی- میلاد صادق یگواری- حسین قائدی 
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قطعات سنگ در مسیر حرکت در طی ریزش انرژی 
جنبشی آن افزایش می‌یابد که بعد از جهش و برخورد با 
سطح شیب مقداری از انرژی جنبشی کاسته می‌شود. این 
کاهش به‌دلیل تغییر محسوس در ضریب ارتجاعی سطح 
جاده نسبت‌به سطوح دیگر شیب در محدوده‌ی جاده 
بیشترین مقدار را دارد. انرژی جنبشی در انتهای مسیر 
به‌دلیل کم شدن شیب دامنه و تغییر نوع حرکت از حالت 
جهشی به غلتشی و لغزشی, با کاهش روبه‌رو می‌شود. 
شکل (۷) انرژی جنبشی در طی ریزش را نشان می‌دهد. 
هم‌چنین مطابق شکل (۸) سرعت انتقالی در طی مسیر 
رفتاری شبیه به رفتار انرژی جنبشی دارد و تغییر نوع 
حرکت باعث نوسانات سرعت می‌شود. 

باتوجه به نمودارهای به‌دست آمده از شبیه‌سازی 
سنگ‌افت می‌توان به این نتیجه رسید که با افزایش وزن 
قطعه سنگ تغییر خاصی در میزان بیشترین مسافت طی 
شده و سرعت. حاصل نشده‌است و انرژی جنبشی کل 
قطعه تنها عاملی است که دست‌خوش تغییرات شده و 
به‌طور آشکاری افزایش یافته‌است. بنابراین می‌توان از 
ویژگی‌های تویوگرافی‌های دامنه و وزن قطعه سنگ‌های 


سال سی و چهارم» شمارة یک» ۱۶۰۰ 


ایجادشده در منشأً افت به‌عنوان عوامل تعیین کننده‌ی نوع 
و میزان متغیرهای مختلف حرکت نام برد. باتوجه به منشاً 
سنگ‌افت و توپوگرافی دامنه‌ی موجود. همه‌ی نقاط میان 
منشاً سنگ‌افت و محل توقف قطعات ریزشی» مستعد 
برخورد قطعات سنگی است که می‌تواند سبب 
خسارت‌های زیادی شود. بدیهی است که برخورد 
قطعات در پهنه‌های بالاتر بیشتر است. ازاین‌رو باتوجه به 
محل توقف مصالح ریزشی می‌توان محدوده‌ی در معرض 
خطر را تعیین کرد. 

مطابق نمودارهای به‌دست‌آمده و محل قرار گرفتن 
جاده و تأسیسات انسانی در دامنه‌ی شیب درصورت 
وقوع سنگافت در منطقه‌ی مورد نظر قطعه سنگ‌های 
جدا شده از منشاً سنگ‌افت با ارتفاع جهش ۷ متره 
شرصق ۱۲ قزر کالید و مین آنرژی ضیف ۱۷۰ 
کیلوژول است و باعث ایجاد خطر ریزش در جاده 
می‌شود. هم‌چنین باتوجه به ماکزیمم مسافت طی‌شده 
سرعت ۲ متر بر ثانیه و انرژی جنبشی ٩۰‏ کیلوژول قطعه 
سنگ باعث به‌وجود آمدن خسارت‌های جانی و مالی در 
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از رسیدن قطعات به جاده. می‌توان از دیوار حائل استفاده 
کرد. باتوجه به نتایج شبیه‌سازی به‌منظور جلوگیری 
سنگ‌افت. مطابق شکل ٩(‏ محل قرار گرفتن مانع 
باتوجه به افت انرژی در بالای جاده در ارتفاع ۱9۷ 


تحلیل بدیده‌ی سنگ‌افت در دیواره‌ی پرشیب بااستفاده از مدل‌سازی... 


جلوگیری کرد. هم‌چنین باتوجه به شکل (۱۰) و ارتفاع 
پرش, ارتفاع مانع حداقل باید ۵ متر درنظر گرفته شود. 
ظرفیت مورد نیاز اين مانع پس از شبیه‌سازی ۲۰۰ 


کیلوژول به‌دست آمل‌هانسست: 


متری می‌باشد تا بتوان از رسیدن سنگ‌ها به مسیر جاده 


۶۴6۲2۷ ۱6۱۳6۲۱6 ا0]" 


سرعت‌اولیه ۱.۲۱ متربرئانیه سس 
یدون سرعت‌اولیه عم 


وزن اولیه ۱۰۰ کیلوگرم 


50 100 150 200 250 300 350 400 


۱۵621100 ]۳[ 


شکل ۷ انرژی جنبشی در طی ریزش (ژول) 


۲۲۵۳05۱8110۴۱۱ ۷۵۱۵۵۱۰۷ ۵ 


سرعت اولیه ۱.۲۱ متربرتانیه سس 


100 150 200 250 300 350 00 
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شکان ۸ سرت اکقال کرک سر (شر نی فانت) 
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[۱۳۰/5] ۷۱۵۵۲۵۲ ادده‌ناهاع۲۳2 


سال سی و چهارم» شماره یک ۱۶۰۰ 


هادی فتاحی- میلاد صادق یگواری- حسین قائدی 


ی اک 
وری غربه‌گبر پرای چاوگیری از حرکت سنگ افت‌هاء 


شکل ٩‏ محل قرار گرفتن مانع 
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شکل ۱۰ انرژی جنبشی کل بعد از اجرای مانع 
نتایج حاصل از شبیه‌سازی به‌روش المان مجزا به‌وسیله‌ی مجموعه‌ای از عناصر کروی‌شکل و باتوجه به 
۳۳3 خواص آنها مطابق شکل (۱۱) تنظیم شده‌است. 


شبه‌سازی مقطع دیواره و بلوک حداشده در تراز با ین 


سال سی و چهارم. شمارة یک, ۱۶۰۰ با میلست عم اف ریس 


تحلیل بدیده‌ی سنگ‌افت در دیواره‌ی پرشیب بااستفاده از مدل‌سازی... 
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شکل ۱۱ شبیه سازی مقطع دیواره و بلوک جداشده در نرم‌افزار ۳۳0۳ 


در شببه‌سازی بااستفاده از روش المان مجزا و 
نرم‌افزار ۳ می‌توان با تحلیل نمودارهای به‌دست‌آمده 
از شبیه‌سازی, به بررسی بیشتر فرایند سنگافت و تحلیل 
فرایند باتوجه به این‌که شبیه‌سازی در سه بعد انجام 
و ذرات با دیواره وجود دارد به‌صورت خودکار شناسایی 
مراحل مختلف را به‌دست آورد؛ به‌عنوان مثال شکل (۱۲) 
در این مقاله. فرایند شبیه‌سازی در ۶۰۰۰۰ مرحله 
برای بررسی حرکت قطعات در طول دامنه انجام 
شندهامست, 

باتوجه به نمودار انرژی جنبشی کل در شکل (۱۳) 
مشاهده می‌شود بلوک پس از جدا شدن از جایگاه ابتدایی 


خود به‌دلیل وزن زیاد بلوک, دارای انرژی جنبشی بالایی 


نشریه مهندسی عمران فردوسی 


تبدیل شدن آن به قطعات کوچک‌تر (شکل ۱۳) به‌شدت 
خافقن ض بالق تفارش از ری خعودفت وه زار 
به‌دلیل برخورد کاهش می‌یابد. پس از این کاهش انرژی 
قطعات حاصل شروع به حرکت در سطح دامنه می‌کنند 
که این حرکت همراه با افزایش انرژی جنبشی است. اما 
این افزایش تا زمانی ادامه می‌پابد که قطعات با برخورد 
به سطح دامنه انرژی و سرعت خود را از دست می‌دهند 
و نوع حرکت آنها از غلتش به غلتش با سرعت کمتر و 
لغزش تغییر می‌کند؛ و يا این‌که قطعات به سطح جاده 
می‌رسند و بعد از برخورد با سطح جاده تعدادی از 
قطعات متوقف می‌شوند. اما بقیه‌ی قطعات به حرکت 
خود ادامه می‌دهند که این حرکت باعت افزایش انرژی 
جنبشی آنها می‌شود ولی این افزایش بعد از مدتی رو به 
کاهش و صفر شدن می‌گذارد و دلیل آن کم شدن شیب 


سال سی و چهارم» شماره یک ۱۶۰۰ 


هادی فتاحی- میلاد صادق یگواری- حسین قائدی 
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شکل ۱۲ محل قرارگیری قطعات بعد از ۱۳۰۰ گام 


5 مدید 
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1-1001( 
۲ ۷61 3 
واعع] 
6+ <> 7.1386+001 


4و ۷5 
معا 
10276005 <2> 1.0006+001 


۰ 6۲۵۷۵۶ وطتااناعصمت فععهاا 
۱۱۸۵ ۱۷۱۱۱ ,عناه۱۷۱۱۸8۵۵ 


شکل ۱۳ انرژی جنبشی کل (ژول) 


انرژی جنبشی دو توپ با محل توقف متفاوت در 
شکل‌های (۱۶) و (۱۵) نشان داده شده‌است. در نمودار 
شکل (۰)۱۶ توپ شماره‌ی ۸ در محدوده‌ی جاده» انرژی 
جنبشی برابر با صفر دارد که نشان‌دهنده‌ی این است که 
در سطح جاده متوقف می‌شود و قبل از آن انرژی جنبشی 
معادل ۶۰۰ کیلوژول دارد؛ نمودار شکل (۱۵) انرژی 
جنبشی توپ ۲۹ را نشان می‌دهد که این توپ با یک سیر 


سال سی و چهارم» شمارة یک» ۱۶۰۰ 


صعودی در انرژی جنبشی حرکت می‌کند و پس از 
برخورد به سطح جاده انرژی آن به‌شدت کاهش می‌یابد 
و با عبور از سطح جاده انرژی جنبشی آن تا رسیدن به 
پایین شیب افزایش می‌یابد و پس از آن انرژی آن کاسته 
می‌شود و در فاصله‌ی ۳۵۰ متری به توقف کامل می‌رسد 
و انرژی جنبشی آن صفر می‌شود. 

همان‌طور که در نمودار شکل‌های (۱1) و (۱۷) 


نشریه مهندسی عمران فردوسی 


مشخص است توپ‌ها (بیانگر قطعات سنگ می‌باشند), 
حالت و مکان توقف‌های متفاوتی را دارند. براساس این 
شکل‌ها مشاهده می‌ شود که قطعات حداشده در مسیر 
دامنه به حرکت خود ادامه می‌دهند و تعدادی از آنها در 
محدوده‌ی جاده توقف دارند؛ هم‌چنین تعدادی از قطعات 
که در جاده متوقف شده بودنده با یرویی که دیگر قطعات 
در حال افت به آنها وارد می‌کنند. به حرکت خود در دامنه 


ادامه می‌دهند و در انتهای دامنه متوقف می‌شوند. علاوه‌بر 


۲060۲2۱۱ 5۳2۲۵0۵0 -عا۲ حول 


تحلیل بدیده‌ی سنگ‌افت در دیواره‌ی پرشیب بااستفاده از مدل‌سازی... 


این. قطعاتی هم بدون اين‌که توقفی در طول مسیر داشته 
باشند به‌سمت پایین دامنه حرکت می‌کنند؛ به‌دلیل عدم 
توقف. انرژی جنبشی در این قطعات افزايش می‌یابد و 
این افزايش باعث می‌شود قطعات مسیر طولانی‌تری را 
نسبت‌به دیگر قطعات طی کنند. به‌علاوه باتوجه به 
شبیه‌سازی انجام‌شده که در شکل (۱0) نشان داده 
نانک شقر نم سانه‌مایی. اتقاق‌افتاده. دز فاصلاق 


افقی ۳ متری از منشا ریزش می‌باشد. 
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شکل ۱۵ انرژی جنبشی توپ ۲۹ 


نشریه مهندسی عمران فردوسی 


سال سی و چهارم» شماره یک ۱۶۰۰ 


هادی فتاحی- میلاد صادق یکواری- حسین قاثدی 
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شکل ۱۷ بیشترین مسافت طی شده توسط توپ‌ها 


براساس نمودارهای به‌دست آمده از شبیه‌سازی 
بااستفاده از نرم‌افزار المان مجزا *۳۳ می‌توان نتایج زیر 
را به‌دست آورد: 
باتوجه به نمودار انرژی جنبشی کل حاصل از سنگ‌افت 
۱ در طول مسیر حرکت مشاهده می‌شود که وقتی 
قطعات به مسیر جاده می‌رسند. انرژی جنبشی آنها 
معادل ۰۰ تا ۱۰۰ کیلوژول است و هم‌چنین قطعات 


سال سی و چهارم» شمارة یک» ۱۶۰۰ 


با انرژی جنبشی معادل ۱۰۰ کیلوژول به مناطق 


خسارت‌های جبران‌ناپذیری در پی خواهد داشت. 


درصورت 


۲. همان‌طور که در نمودار انرژی جنبشی مشخص است. 
وقتی بلوک به‌صورت یکپارچه شروع به حرکت 
می‌کند انرژی جنبشی زیادی دارد و وقتی به قطعات 
کوچک‌تر تبدیل می‌شود انرژی آن به‌شدت کاهش 


نشریه مهندسی عمران فردوسی 


۱۱۲ 


یز ان ان ارو نی من 
رسید که کاهش ابعاد و وزن قطعات یک روش مناسب 
برای کم کردن انرژی جنبشی قطعات» خسارات و 
خطرات ناشی از آنها است. 

۳ باتوجه به نمودارها و شکل توقف توپ‌ها در دامنه. 
اگر سنگ‌افت اتفاق بیافتد. تعداد زیادی از قطعات در 
محدوده‌ی جاده‌ها و مناطق مسکونی توقف دارند که 
این توقف‌ها راه‌بندان در جاده و خطرات مالی و جانی 
را در پی دارند. 

6 باتوجه به این که بیشترین فاصله‌ی افقی طی‌شده توسط 
قطعانت ۱۳۹۲۰ رم فرتناشان: ما معا وت ور ان 
سنگ‌افت بسیار بالا می‌رود چرا که فاصله‌ی افقی 
مناطق مسکونی ۲۶۰ متری می‌باشد و تعداد زیادی از 
قطعات به این فاصله خواهند رسید. 

برای جلوگیری از رسیدن قطعات به سطح جاده 
همان‌طور که گفته شد یکی از روش‌ها می‌تواند کم کردن 
وزن و ابعاد قطعات در منشاً ریزش باشد. روش دیگری 
که می‌شود از آن استفاده کرد ایجاد یک دیواره‌ی حائل 
در مسیر است که باتوجه به نتایج حاصل از نرم‌افزار 

0/1 می‌توان محل قرار گرفتن آن را مشخص کرد. از 

طرفی چون قطعات در مسیر حرکتی خود به سطح جاده 

برخورد می‌کنند و بعد از این برخورد به‌صورت جهش به 
حرکت ادامه می‌دهد. برای جلوگیری از ادامه‌ی این 
حرکت می‌توان از جوی ضربه‌گیر که کف آن با ماسه 
پوشیده شده و یک دیواره در کنار جاده استفاده کرد. اين 
کار باعث می‌شود انرژی قطعات در برخورد با ماسه 
کاسته شود و بااستفاده از دیواره‌های فلزی یا بتنی کنار 
جوی از ادامه‌ی حرکت آنها جلوگیری کرد و قطعات را 
در کنار جاده به دام انداخت. باتوجه به این که توپوگرافی 
شیب عامل مهمی در سنگ‌افت به‌شمار می‌رود. تغییر در 
توپوگرافی با ایجاد یک پله در دامنه نیز می‌تواند روش 
کارآمدی برای جلوگیری از سنگ‌افت باشد. 


نتیجه گیری 


نتایج به‌دست‌آمده از این تحقیق را می‌توان در موارد زیر 


نشریه مهندسی عمران فردوسی 


تحلیل بدیده‌ی سنگ‌افت در دیواره‌ی پرشیب بااستفاده از مدل‌سازی... 


خلاصه کرد: 

. از جمله عواملی که در ارزیابی پدیده‌ی سنگ‌افت تأثیر 
دارند می‌توان به عوامل توپوگرافی نظیر هندسه‌ی 
شیب وزن قطعات سنگی, منشأً ریزش سنگ‌ها و 
خواص مکانیک سنگی دامنه اشاره کرد. باتوجه به 
نتایج به‌دست آمده می‌توان گفت اصلی‌ترین عامل 
کنترل کننده‌ی حرکت قطعه سنگ برروی شیب. عامل 
توپوگرافی است. 

۲ باتوجه به این که انرژی قطعات سنگی با میزان مخرب 
و خطرناک بودن آن رابطه‌ی مستقیم دارد. از این‌رو 
هر چه انرژی جنبشی در قطعات بیشتر باشد میزان 
خسارات ناشی از آن بیشتر خواهد شد. بابراین 
باتوجه به نتایج به‌دست آمده در شبیه‌سازی‌هاه می‌توان 
با ایجاد تغییرات مناسب در عواملی که در سنگ‌افت 
مثر هستند میزان خسارت‌های احتمالی از آن را 
کاهش داد يا به صفر رساند. ایجاد تغییر در هندسه‌ی 
شیب. تغییر در ابعاد قطعات. ایجاد موانع در مکان 
مناسب در مسیر حرکت قطعات از راهکارهای مناسب 
به‌شمار می‌روند. 

۳ در این مقاله سنگافت بااستفاده از نرم‌افزار 00811 
شبیه‌سازی شده‌است و باتوجه به نتایج حاصل از آن 
این پدیده مورد بررسی قرار گرفته‌است. از آنجایی که 
این نرم‌افزار فقط مسیر حرکت را شبیه‌سازی می‌کند و 
عواملی نظیر شکل. اندازه و وزن قطعات در 
شبیه‌سازی تأثیری ندارد و فعل و انفعال بین ذرات در 
مسیر حرکت درنظر گرفته نشده‌است و ذرات را جدا 
از هم و منفرد در نظر گرفته‌است و باتوجه به این‌که 
در واقعیت قطعات در حال ریزش و سقوط برروی 
یکدیگر تأثیر گذارند» نتایج حاصل از شبیه‌سازی این 
نرم‌افزار با واقعیت سنگ‌افت فاصله دارد. از این‌رو 
برای پوشش و رفع نواقص و ضعف‌های این نرم‌افزار 
از نرم‌افزار ۳۳8 استفاده شده‌است. در این نرم‌افزار 
علاوه‌بر در نظر گرفتن فعل و انفعال بین ذرات و 
دیواره‌هاء فعل و انفعال بین یک ذره با ذرات دیگر نیز 


سال سی و چهارم. شمارة یک» ۱۶۰۰ 


هادی فتاحی- میلاد صادق یگواری- حسین قائدی 


در نظر گرفته شده‌است تا به شرایط موجود در طبیعت 
نزدیک‌تر باشد تا بتوان نتایج قابل قبول و مطمئن‌تری 
را از شبیه‌سازی به‌دست آورد. 

. بااستفاده از اين دو نرم‌افزار که مکمل هم هستند 
می‌توان تحلیل بهتری از سنگ‌افت در یک دیواره 
داشت و براساس نتایج حاصل از آنها خطرات 
احتمالی در صورت سنگافت را مشخص کرد و 


۱۳ 


درجهت کاهش این خطرات انجام داد و خسارت‌های 
جانی و مالی ناشی از آنها را به حداقل رساند. 


. ضمنا پيشنهاد می‌شود به‌عنوان تحقیق بعدی برای 


اندرکنش داخلی توده در نرم‌افزار 0 مجموعه‌ای از 
دانه‌ها (با اندازه‌ی کوچک‌تر) و انتخاب پارامترهای 


واقعی‌تری از سفوظ سنگ دو:مطفقه به‌دست. آیت, 
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1 ۶ وفطممت0م۷ افممتمهام۴! 0ص ممتبمتاعمک ۵1 اممه‌تاهم رل پ۷۷۵ 200 ویک ر۷۷۵۵8۵ ,.؟ .ک بتالفط) 
م۱ ۱۱۰۰۱۰۵۱۱۵ 
2002(۰) 

40 م ۱۷۲020 2 عمماظ عملمن۱۷]۵ ۳۱۵2۵۲0 ۴5۵۵1۵11 ,۷۲ بلتعماته9 220 ب ,1200760011 .۲ ,۳11111270 
۴1 ۵۲۱۱ ۵۱۱4 و۵20 ۸۷۵۸۵ , ۸۵۲02 عصتافک؟ تمامحصصهت2ظ 19-۵260 2 عصتونا 9۳/1۱267240 ظ1 
2003(۰), 435-442 .00 5 0۰ ,3 ,۷۵۱ ,96011665 

2 ۸۱۵92 ا]ممصومع‌عمه علفنک هه ۳22274 ۵121م۴. ریق بتعتات) 20 ررظ بطهطصمطمامط رب متااع00227 
۰ 34 .۲۷۵1 ,۱۱۵/۵8۵6160 ۲۱۷۱۳۵۵۱۱۵۵۵ ۱۱21۷ تحعاصمی رماله ۷ ۱۲۵۲۵ عظ) ا عم من صمم‌وهصهنداز0 
2004(۰), 191-208 .00 ,2 

۵08 ۳08۵ 10۲ ۱2۵ ۳۱۵2۵۲0 ۴۵۵۲211 60ع2ظ وان وا بک ,۷/۵0۵ 200 .لا معظ2؟ و.؟ ک بتافط) 
٩۱7۳۳, 1 ,)2004(۰‏ .ار 41 ۷۵۱۰ رکصع02 ۱/۱۱۱۵ ۵ وع۷۵۵۱:۵( ۵۵۵1 ال ۲۱۱۱۵۳۱۱۵۵0۵ 

صما )وع:۲۵0 2 طا 0۷۱0اع۸ هم امومع مه لهتادمو وبط ر2امطمعنک م2 رب ,الما .و ,۳۵۳۲ 
۰ ,1-4 .۱۱0 ,74 ۷۵ وما۵0۵0۵۵ع رتاو عفمت ل2من۵102ظ1(60۲02601۴80۲۵ م--و۱۵۵۱۵ظ ۳/195 معط 
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۱۴ تحلیل بدیده‌ی سنگ‌افت در دیواره‌ی پرشیب بااستفاده از مدل‌سازی... 


219-231 ,)2006( 

00 ۲اتزاتانمعویاو متباته۳ ۵۶ احمصصووعوفم. .نا ۷۲۷۰ ,متصاظ مضه رب بوع۵۱۷ 016112 42 ,.۲, 17 520072 
2006(۰), 211-224 .00 4 .۱0 ,86 ۷۵۱ جوهاههع و9۱۵2 ,۱0۵۵16 ۲۵22۷ عصلونا 510065 

8 ۳۵2۵70 0601021021 0۶ ون مطلووعووه مه عصلم۱۷]2۵. ر.نا رتم12۱ 0صه .و رعلقتت۷]۱۵ ,۷ بطموتعاظ 
۵5۵061 6۱۵۲۵۲۱۲ ,"نامع رتمالج ۷ بعطامصت۲ مج هابافمتصعظ مملوخ )مععطها عط 
۰ ۱۱0۰ رع۵00ک ,5۱/۳۷6 ۲60108661) ,6509۲6۵6 0۴۸۷۵۱۵۵۵ ۵۵۵۲۱۱۵( را وت 0۷۵/۵۵۵۱۵۵۱۵ 
2006(۰), 06-1 

16۷ ,029016 ۸0 0صنامتج ۲۱۵2۵۲0 ۴۵۵1211 ۵۶ اممجطووعوو۸ ۰ به .نا مع020 220 ,۷ بصتاظ و.] م100۵ 
.(2007), 191-200 .00 و1 0۰ ,53 ۷۵۱ ر(ع0/0ع6 ۳۱۱۷۱۲۵0۳۵۱۱۵۱۱۵ 

0ص لمممزممممونل-ععنط ]۰ 1۵۲ ممتقصم ۱ ورن ۸ :)دمص هم بلط ) بصطشا 224 وناب بطلات۷2 بط مطش] 
-262 .00 و2 ٩/0.‏ ,33 ۷۵ 660536161665 ع 0۳۱۵۵/۵۲5 ,۱۷۲۵061182 ۳۱۵۸2۵۲0 ۹۵۵1۵11 0ماباطازوزا رللهتاهم5 
2007(۰), 279 

۷/1۴ ]1۳۵۵0 ج0 فلله۲ مک ۵۶ مماه)صمصصع ۳۳۵ ده فمتلبای لقاممصهص ۳ ریله اه ریق رتصتطمعهزن 
و 708-715 .00 ,4 ۱0۰ ,46 .۵1 ۷ ,و66 ۷۵۱۵ ۵۳0 ۵۴۵۵1۷۵6۱۵066 ۵۱۱۵ ۲۱۱۵۳۱۱۵۱۵۱۱۵ 5۱۲۵۵69 
2009(۰) 

عط ۶ تاد ععم ۸ :ههام۲۶ عم ماما ممتاموناوه 1 آم‌نومامعن0) مصلرممصامص. یله ۵ .لا مدفط) 
,150۲6۵ 1۳ مر ۳۱۵۲1۲۵۵6 ۷۷۵۲۱۵ وتا معط اه ۲2۵۵ عم ۵ ۵۵ 02۳۷۵0 معفحط ۹601822760۶۲۵6 5962160 
2009(۰), 69-78 00۰ و ٩۵۰‏ و13 ,۷۵۱ رال 6605661665 

هام۲ امدآ مج هلت 2 طظا کممصصووعوفه ۳۱۵2۵270 ۴5۵۵1911 ریک ب1ع20 224 ,۳ بتالعمنالوناهن ]1 
2009(۰), 963-972 00۰ و5 ٩0۰‏ ,56 ۷۵۱۰ ,60108 ۲۱۷۱۲۵۳۸۵۱۱۱۵۱۸۵ ۲۱۵۷(۲۱ رهله0200200) ر029016 تهولط0]2) 
,1512040 تصعمنصهی ب۲۵۲ ومتصتطته 10۲ مه للااتامویوته ال2گهم۴. ر.ط رففااما 2020 بظ ظ بتامتطماحظ 
٩۲0. 1-2, ۵0. 152-166 ,)2010(‏ ,118 ۷۵ ,0660۳71070۳00 ,01۲66686۲ 

۵ ,۱۷۲0۵061۲ ظ۵۵ظ ]۳۵۲۵۹ 2 صنط۳/1 ۱۷۱۵0016 ۴۵۵1211 ب3-1 2 ۵۶ ممتامناله۳۲۷. مبلق 6 و۱۷۷۰ رتعحصصصره ۴[ 
2010(۰), 699-711 .00 ب4 .۱0 ,10 ۷۵۰ ,۹6۵665 6۹/۵۱ ۵۳۱۱ ۵۱0 ۳۵20۲05 

6 0۵۱10۷ 10۲ متعایتت ااحمصتقک که ممیهعنامع۲ ۱ اماممهن ۳ رن رصتد ۷ 20 نا رقصهتا]۲ 
(2010), 326-338 .00 ,3-4 .0 ,120 ۷۵۰ ,چع660۱10۳۳0/۵0 , وعتاآن۳۵ 

ک 20 ۸۷۵۸۸۵ ۹۵۷۱۵۷۷۲۱ ۸ -ط10ا0ع۳۳۵]6 هام5 مه مماه2120)۵18) ۰۴۵01211 ریلح اه ربظر طله01167 ۲۷ 
.(2011), 2617-2651 .00 ,9 ۱0۰ ,11 ,۷۵۱ ,6۵665 5۹/۵ ۲۲۱۲ ۵۵ 

فعصتاهگ 0601021621 1۳11۴6۲6۵۲ 1۵۲ ومعصه۴ ۱6۲2۷ ااهزمم۴ ۵۶ ممتاماه سا لهمتادتاهاو یله اه رب رتته5020 
٩0. 00. 57-65 ,)2013(‏ ,158 ۷۵۱۰ ,60و ۲۱۱8/۱۱66۲۸8 ,"2112داخ روعه ۷۷ طانامو سعل ۶ 

01 0۷ طملله۱ رمتقمتما۱ تومالخ 26 22270 ۴۵6121 ۵۶ اممصووععف۸. یله ۵ با ۸ بلععفتا0 ۷ 
۴0 0۴ 9۲۱۱۵ هل ,۱۵۷۵۹0۵۵0108 ۳۱۵۱0 هه دنا مدع عمجم لهع۵مصه-مهم عصلونا (مصادتض 
(2015), 309-321 .00 .0 ,101 ۷۵۱۰ ,9667۱6 ۵۲۱/۱ 

و۱0۵۷ ۵207۵۷ وتات 02اناهمم2 ۵۶ لهتامعاهظ الفتهمع ۶ ممتاهعتاوه۳۷]. وله له ر.ظ مطلادمکز 
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11۰ 


12 


13. 


14 


13: 


16۰ 


17 


18. 


19۰ 


20 


21۰ 


22. 


23. 


24. 


هادی فتاحی- میلاد صادق یگواری- حسین قائدی ۱۱۵ 


۸۲۵۵ 0۲۵ 0۴ 66056۵0665 ۷۵۱ 13, ۱0. 16, 00. 1-11 ,)2020( 

۱160انددرمدط فصته)صیم۱۷ 0عم)مطو ما )م۱ متصصولهو افو گ0 فلو افص ۱۷۲۵۵6 فتلته۳ یله ام ریک رفظ .25 
.(2021), 425-446 .00 م1 .۱۱0 ,18 ۷۵۱۰ ,۲۵۱05/105 ,۱۷۲۵0611182 آهعزرمجرب آ 0ج متامع ]18۷65 16۱۵1160 

0 مومع امعم ۵۶ معلوعنز رولوزلفممه ملافتازطاهطا۳:۵ ۵ 100 ۸ بالع۴۵0۲۷ .نا ۷۷۰ ,قطع۵7٩‏ .26 
1998(۰) ,۲۴۵612115 0۶ ۳۲۱601610 

:7 ۳۳۵/۵۲۱۱ 0006 ۲۱۵۷۷ ۳۵۲۵۵۵ عط) عملونا وال۲2 وفه۱۲ ۴۵۵ ۵۶ وصتاام0هی ناد ری بله‌عژ۳۲0 220 ۸ بط۵عظ۳ .27 
.102 

عصلوت] ممتاهع6۳1 ۷ هه 5126 صنقری مه عصلل‌صهم1(6 کون معصمل‌نوهان ٩‏ مصضا 0۶ ممتامتااه ۲ ریله اه .1 ,عم .28 


2017(۰), 133-141 .00 ,2 .0 ,27 ۷۵۱۰ و20 0۴1۳۱۵۱۵۵۳۱8۵ 0۵۱۵۲۵ 16 ,ولو رادم لهعزمصیل 


سال سی و چهارم» شمارة یک, ۱۶۰۰ نشریه مهندسی عمرال فردوسی 


۱۶ تحلیل بدیده‌ی سنگ‌افت در دیواره‌ی پرشیب بااستفاده از مدل‌سازی... 


نشریه مهندسی عمران فردوسی سال سی و چهارم» شماره یک ۱۶۰۰ 


